Small regulatory RNAs may sharpen spatial expression patterns 
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The precise layout of gene expression patterns is a crucial step in development. Formation of a sharp boundary 
between high and low expression domains requires a genetic mechanism which is both sensitive and robust to 
fluctuations, a demand that may not be easily achieved by morphogens alone. Recently it has been demonstrated 
that small RNAs (and, in particular, microRNAs) play many roles in embryonic development. While some RNAs 
are essential for embryogenesis, others are limited to fine-tuning a predetermined gene expression pattern. Here 
we explore the possibility that small RNAs participate in sharpening a gene expression profile that was crudely 
established by a morphogen. To this end we study a model where small RNAs interact with a target gene and 
diffusively move from cell to cell. Though diffusion generally smears spatial expression patterns, we find that 
intercellular mobility of small RNAs is actually critical in sharpening the interface between target expression 
domains in a robust manner. We discuss the applicability of our results, as examples, to the case of leaf polarity 
establishment in maize and Hox patterning in the early Drosophila embryo. Our findings point out the functional 
significance of some mechanistic properties, such as mobility of small RNAs and the irreversibility of their 
interactions. These properties are yet to be established directly for most classes of small RNAs. An indirect yet 
simple experimental test of the proposed mechanism is suggested in some detail. 



INTRODUCTION 

Morphogenesis proceeds by sequential divisions of a de- 
veloping embryo into domains, each expressing a distinct set 
of genes. Each combination of genes is associated with a 
particular cell identity. At advanced stages of development, 
most genes that define cell identity are either highly expressed 
('on') or strongly inhibited ('off') in a given cell. For ex- 
ample, two adjacent domains may be differentiated by a high 
expression of some genes in one, and low expression in the 
other. In such cases it is important that cells of the two pop- 
ulations do not intermix. Moreover, the number of cells that 
show intermediate levels of expression, typically found at the 
interface between the two sets, should be kept to a minimum. 
These demands are necessary in order to unambiguously de- 
fine the identity of each cell. A spatial gene expression pattern 
that obeys these demands is said to exhibit a sharp interface. 

A crucial step in setting the interfaces of gene expression 
patterns is often the establishment of a concentration gra- 
dient of molecules called morphogens. Some morphogens 
are transcription factors which regulate gene expression di- 
rectly fHIIl]- Others are ligands that bind cell-surface recep- 
tors signaling the activation of target expression j^]. Since 
morphogens act in a concentration-dependent manner, a mor- 
phogen gradient may be transformed into a gradient of its tar- 
get messenger RNA (mRNA). 

In principle, a single morphogen interacting cooperatively 
with its target enhancer can generate a sharp interface in 
the target transcription profile, by modulating the rate of its 
mRNA transcription as a function of the nuclear spatial coor- 
dinate 0j, This may be done e.g. by cooperative binding to a 
receptor or to a promoter |5[ or by zero-order ultrasensitivity 
As an example, in fruit-fly early embryonic development 



Hunchback transcription depends on the cooperative binding 
of about 5 Bicoid molecules 0]. An obvious limitation in this 
mechanism is the need for large cooperativity factors or cas- 
cades of reactions, which make it prone to fluctuations and 
slow to adapt ansa. Recently, the role of small regula- 
tory RNAs in establishing developmental patterning has been 
observed in plants \VA [Till and animals JTill . In particular, 
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it has been suggested that microRNAs (miRNAs) confer ac- 
curacy to developmental gene expression programs 11511 . This 
raises the possibility that small RNAs aid morphogen gradi- 
ents in establishing sharp interfaces between 'on' and 'off 
target-gene expression. 

In this study we formulate a mathematical model in which 
small regulatory RNA help morphogens to determine cell 
identity by sharpening morphogen-induced expression pat- 
terns. For specificity we assume here that the small RNA 
belongs to the miRNA family, and consider other classes of 
small RNA in the Discussion. miRNAs constitute a major 
class of gene regulators that silence their targets by binding 
to target mRNAs. In metazoans primary miRNA transcripts 
are transcribed and then processed both inside and outside 
of the nucleus to form mature transcripts approximately 21 
nt in length that are then loaded into the RNA-induced si- 
lencing complex (RISC) 111 611 . They are found in plants IU7I1 
and animals 1 1_8], including human lfl9ll . and are predicted 
to target a large fraction of all animal protein-coding genes 
1 19[ Ufx Gill . In plants miRNAs are known to affect mor- 
phology 111 U Il2l 12211 implying that they play an important 
role in determining cell identity. This is underscored by the 
fact that the spatiotemporal accumulation of miRNAs is under 
tight control in plants B23I1 . fly B 14112411 and zebrafish 1250 . 

Our model is constructed in one spatial dimension, namely 
along one spatial axis. Domains of gene expression are laid 
out along this axis, and we assume no significant variance 
along other, perpendicular axes. Two key ingredients of the 
model are a strong interaction between miRNA and mRNA 
and intercellular mobility of the miRNA. Within this frame- 
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work miRNAs generate a sharp interface between those cells 
expressing high levels of the target mRNA and those express- 
ing negligible levels of mRNA. We use physical arguments 
to understand the range of parameters where this sharpening 
occurs. Our model predicts that the spatial position of the in- 
terface is precisely determined: mobile miRNAs spatially av- 
erage individual cellular fluctuations without compromising 
the interface sharpness. We use computer simulations to show 
that this is also true even with low copy numbers. A conse- 
quence of our model is that a local change to the transcription 
profiles can induce a nonlocal effect on the mRNA concentra- 
tion profile; we outline an experiment to detect this nonlocal 
property. We also consider the special case where the miRNA 
and target transcription profiles coincide lEill . We show that 
miRNAs preserve their sharpening function provided mobile 
miRNAs change irreversibly to an immobile state before in- 
teracting with their targets. By separating mobility and inter- 
action in this way, the target is rescued from miRNA-mediated 
silencing in selected cells thereby generating a sharp expres- 
sion pattern. Finally we consider possible applications of 
these ideas in plants and fruit fly. 



MODEL AND RESULTS 
Formulation of model from available biological evidence 

Our theory comprises three central elements. First, miRNA 
and its target are taken to be transcribed in a space-dependent 
manner. Second, we assume that the interaction between 
miRNA and target irreversibly disable the target mRNA from 
being translated into proteins; this for example may be done 
by promoting the degradation of the target. Furthermore, the 
miRNA molecule itself may be consumed during this interac- 
tion. Last, we allow for the possibility that miRNAs move 
between cells. Before defining the model, let us review the 
available data regarding each of these processes. 

microRNAs and their targets are often expressed in a coor- 
dinated manner l27ll . We consider two classes. In the first, the 
regulatory network is designed to express the miRNA and its 
targets in a mutually exclusive fashion. For example, the ex- 
pression patterns of the miRNA miR-196 and its target Hoxb8 
are largely non-overlapping in mouse II26I1 and chick 12811 . 
Similarly, the nascent transcripts of Ubx, still attached to the 
DNA, are expressed in a stripe near the center of the early 
embryo, while nascent transcripts of its regulator, iab-4, are 
simultaneously observed in nuclei posterior to this domain 
1 2911 . A recent large-scale study in fly showed that miRNAs 
and their target genes are preferentially expressed in neigh- 
boring tissues 111 511 . Likewise, in mouse 113011 and in human 
1 3 If! predicted miRNA targets were found at lower levels in 
tissue expressing the cognate miRNA than in other tissues. 
Cases in the second class are characterized by co- 



gene cluster may target co-expressed Hox genes (26 

Our model assumes that the synthesis rate of the miRNA 
and its target are smoothly varying along a spatial axis, x. 
This, for example, may be the result of a common morphogen 
regulating (either directly or indirectly) the two species. The 
transcription profiles a^(x) and a m (x) of the miRNA and its 
target, resp., may be anti-correlated (the first class above) or 
highly correlated (the second class). 

The detailed interaction between miRNAs and their targets 
is currently a topic of intense investigation J34, 35]. miRNAs 
induce the formation of a ribonucleoprotein complex (RISC). 
Targeting of a specific mRNA by a RISC is done via (often im- 
perfect) base-pair complementarity to the miRNA j 18J]. Upon 
binding, protein synthesis is suppressed by either translation 
inhibition or mRNA destabilization 134LI35I1 . While it is likely 
that miRNA can go through a few cycles of mRNA binding 
1 3611 . the increased endonucleolytic activity conferred by the 
miRNA makes it plausible that the miRNA is sometimes de- 
graded in the process. In addition, evidence suggests that 
mRNAs which are translationally repressed by miRNA may 
be co-localized to cytoplasmic foci such as processing-bodies 
1 34[ l35l 137 L which are enriched with endonucleases, or stress 
granules 13811 . Taken together, these facts make it improbable 
that miRNAs act in a fully catalytic manner. 

A pair of mRNA-miRNA reactions that describe a spectrum 
of plausible scenarios is 
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where m represents the mRNA concentration and /1 repre- 
sents that of the miRNA. Here, 8 is the average number of tar- 
gets degraded by a given miRNA before it is itself lost in the 
process. These reactions may be realized in different ways. 
For example, the two species may reversibly form a complex 
that is then subject to degradation. Another possibility is that 
the two species irreversibly associate to form an inert com- 
plex. Furthermore, the reaction between the species may be 
reversible, as long as the typical dissociation time is much 
longer than the relevant biological time scale. One way in 
which the cell may control the dissociation time is by regulat- 
ing exit of the RNA pairs from P-bodies D9I1 . 

Can miRNAs move from cell to cell? siRNAs, another im- 
portant class of small RNAs, are known to illicit non cell au- 
tonomous RNA silencing in plants, worms, fly and possibly 
mouse (reviewed in J40I1 ). This may also be the case for trans- 
acting siRNA 111 311 . Evidence in favor of intercellular mobility 
of miRNA is found in pumpkin 114 111 . There, miRNAs have 
been found in the phloem sap which is transported throughout 
the plant by phloem tissue. In animals, many small RNAs, in- 
cluding many miRNAs, were found in exosomes from mouse 
and human mast cell lines, which can be delivered between 



transcription of the miRNA and its targets 12711 . as occurs for 
example in mouse cardiogenesis 0211 and human cell prolif- 
eration 13311 . Likewise, microRNAs encoded within the Hox 



cells J42]. 

In our model we consider the possibility that miRNA mi- 
grate from cell to cell. Mobility of the miRNA species is 
likely to rely on active export from the cell followed by im- 
port to neighboring cells, or perhaps on direct transport be- 
tween cells. On the tissue scale, these transport processes are 
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expected to result in effective diffusion. We therefore ignore 
the small-scale transport processes and model miRNA mobil- 
ity as pure diffusion. 

Finally, we combine these processes into a steady-state 
mean-field model given by 



= — (3^p — kmp, + Dp," . 



(2a) 



(2b) 



The j3 terms describe autonomous degradation (i.e., by pro- 
cesses independent of the other RNA species) and the k term 
describes coupled degradation of both RNA species. Note 
that the case 8 ^ (where miRNA may go through multiple 
rounds of interactions with target mRNAs) can also be brought 
into this form by rescaling ( f2ab . Mobility of the miRNA is 
described by an effective diffusion constant D. The spatial 
coordinate x measures distance along one dimension of a tis- 
sue. All our numerical results shall be presented in units of 
the tissue size, i.e. < x < 1. 



microRNAs may sharpen domain boundaries. 

As described above, a desired target protein profile com- 
prises a domain of cells which express this protein abundantly, 
adjacent to a domain of cells where this protein does not ac- 
cumulate. Furthermore, one requires that the number of cells 
with intermediate expression levels lying in between the two 
domains be minimized — this is our definition of a sharp inter- 
face. In this section we discuss one scenario where the mutual 
consumption of a diffusive miRNA and its target leads to such 
a sharp interface in the mRNA profile. 

We assume that some morphogen controls the transcription 
rate of the target. The transcription profile — namely, the tran- 
scription rate as a function of the spatial coordinate of the 
nucleus — is laid down as a smooth gradient, falling from one 
end of the developing tissue (which, for convenience, will be 
called 'left') towards the other end ('right'). Motivated by 
a recent study that showed that miRNA and their targets are 
preferentially expressed in neighboring tissues 11511 . we focus 
on the scenario where the miRNA transcription is controlled 
in a fashion opposite to that of the mRNA: miRNA transcrip- 
tion is peaked at one end of the tissue (where the mRNA tran- 
scription rate is minimal), and decreases towards the other end 
(where the mRNA transcription rate is maximal). The kind of 
'mutually exclusive' transcription profiles we have in mind is 
depicted in Fig. UK- In this figure, and hereafter, we denote the 
mRNA transcription profile by a m (x) and the miRNA tran- 
scription profile by a^(x), explicitly noting their dependence 
on the spatial coordinate x. We note that, in the absence of 
miRNA-target interaction and of miRNA diffusion, the con- 
centration profiles of mRNA and miRNA simply follow their 
transcription profiles, Fig. [TJ\ and D. Each is rather smeared 
and overlaps the other near the centre of the tissue. 



Before studying the full model, it is instructive to consider 
first the interacting system in the absence of miRNA diffusion. 
In the context of mutually exclusive transcription profiles, we 
expect that each cell would be dominated by one RNA species 
(either the miRNA or the target mRNA), which we will call 
the majority species, and be depleted of the other, minority 
species. In other words, we are making the critical assump- 
tion that the decay of the minority species in each cell is gov- 
erned by the interaction with the other species (rather than by 
its autonomous degradation). This assumption can be made 
quantitative in terms of the model parameters; see Eq. IS4l of 
the Supporting Information. Under this assumption, which we 
refer to as the strong interaction limit, it is straightforward to 
show (Eq.[S6]of the Supporting Information) that the density 
of the majority species in each cell is proportional to the differ- 
ence between the two transcription rates in that cell, whereas 
the minority species is essentially absent. Consequently, in 
the context of mutually exclusive transcription profiles, the 
mRNA level becomes vanishingly small in any cell for which 
Oi m {x) < a^(x), namely every cell to the right of the point 
where the two transcription profiles are equal. The concentra- 
tion profiles of the two RNA species are shown in Figs. [Tj3 
and E, where one can see that the mRNA and miRNA spatial 
expression domains are now complementary and more sharply 
defined. 

The threshold response that arises when both RNA species 
do not diffuse from cell to cell provides insurance against 
the possibility that the mRNA transcription profile is not as 
step-like as is required for unambiguous cell differentiation. 
In other words, miRNA regulation acts as a failsafe mech- 
anism whereby incorrectly transcribed low-abundance tran- 
scripts in the region a m (x) < a^x) are silenced, while 
correctly transcribed high-abundance transcripts in the re- 
gion a m (x) > a^(x) are only mildly affected Il5l 12711 . 
This threshold response in the target profile has been ob- 
served in the context of small RNAs — another class of post- 
transcriptional regulators — in bacteria (43 j. 

We now return to our full model, which allows for diffu- 
sion of the miRNA. To simplify the analysis let us keep the 
strong-interaction limit, described above (and in Eq.[S4]of the 
Supporting Information). In general, one expects that diffu- 
sion makes the miRNA profile more homogeneous, and this is 
confirmed by exact numerical solution of the model, as shown 
in Fig. QJ1 Surprisingly however, the mRNA profile does not 
become smoother. In fact, Figs, [if and F show that this pro- 
file actually develops a sharper drop from high to low mRNA 
levels than there was in the absence of diffusion. More specifi- 
cally, miRNA diffusion creates an interface between high and 
negligible target expression. Increasing diffusion moves the 
interface deeper into the mRNA-rich region and thereby ac- 
centuates the drop in mRNA level across the interface. Al- 
though some miRNA diffusion is required to establish a sharp 
interface in the mRNA profile, the diffusion constant cannot 
be too large. As Fig. QJj demonstrates, increasing the diffu- 
sion constant may result in smearing the interface. A corre- 
sponding increase in the interaction strength, k, can compen- 
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FIG. 1: Sharpening the target expression pattern. A. Transcription profiles of a miRNA (red) and its target (green). Functional expressions 
and parameter values are given in the Supporting Information, unless otherwise noted. B. Steady state concentration of target-mRNA (green) 
and cognate miRNA (red). Regulation by immobile miRNAs (where the diffusion constant is D — 0) removes residual mRNAs from the right 
domain, creating a sharp boundary near the center of the tissue. C. Mobility of miRNA {D = 0.01) further sharpens this boundary by inducing 
an interface between domains of gene expression. D-F. Steady-state concentration of mRNA (green level) for a two-dimensional generalization 
of the model, where diffusion occurs equally in both directions and transcription rates do not vary along the vertical axis. Parameter values 
correspond, respectively, to panels A-C. G. Strong miRNA diffusion (D increased 100-fold) smears the interface between the domains, but 
does not affect the interface position. H. The interface structure is unaffected by a increasing both the miRNA-mRNA interaction k and the 
diffusion constant D 100-fold. 



sate for the increased diffusion and regain the interface sharp- 
ness (Fig.QjT). We will quantify these observations below. 

Diffusing miRNAs can find themselves in one of two very 
different regions. In the miRNA-rich region (including the re- 
gion to the right of the point where the transcription profiles 
are equal) miRNA decay mainly via processes independent 
of their interaction with the target. In this region our model 
boils down to a simple diffusion process accompanied by lin- 
ear decay. Such processes are characterized by a length scale, 
denoted by A, which essentially measures how far a miRNA 
can travel (due to diffusion) before being consumed (by au- 
tonomous degradation). It is thus an increasing function of 
the diffusion constant D, but a decreasing function of the in- 
dependent decay rate /3 M . On the other hand, in the mRNA- 
rich region, a miRNA decays mainly via co-degradation with 
its target. In this region miRNAs decay faster, and one expects 
them to be able to travel over much shorter distances than in 
the miRNA-rich region. In fact, diffusion in this region is 
characterized by another, smaller, length scale, denoted by t, 
which again increases with D but is now a decreasing function 
of the interaction strength, k. Explicit expressions for the two 
length scales are given in the Supporting Information (Eqs.lS8l 
and[S9]>. 

To obtain a sharp interface in the mRNA profile, miRNAs 



should be able to travel from the miRNA-rich zone into the 
mRNA-rich zone. This means that the first length-scale, A, 
should be of the same order as the tissue length. This, for 
example, can be achieved if the diffusion constant D is large 
enough. On the other hand, the vicinity of the interface is 
governed by the other length scale t. This length scale is what 
determines the 'width' of the interface, namely the number of 
cells which exhibit intermediate levels of mRNA expression 
(see blue box in Fig. A sharp interface therefore means 
a small value of I and one way to achieve a small value of I 
is to make the diffusion constant D small enough. These two 
contradicting requirements on D suggest that there might be 
an intermediate range of values for the diffusion constant that 
allows for a sharp interface, but also raises the suspicion that 
this range may be very small and require some fine tuning. 
This , however, is not the case: the fact that A is strongly 
dependent on (3^ (while I does not depend on (3^ at all), and 
that i strongly depends on k (while A does not) means that the 
range of allowed values of D can be set as large as needed. 

In the Supporting Information we develop an approximate 
analytical expression for the mRNA profile in terms of the 
various parameters and the 'input' profiles a m (x) and a^{x). 
There are two lessons to be learned from this exercise. First, 
the interface established by the mRNA-miRNA interaction is 
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effectively impermeable to miRNA diffusion in the strong- 
interaction limit. The system thus separates into two parts 
which — in steady-state — do not exchange molecules between 
them. This property allows one to calculate the position of the 
sharp interface in the mRNA profile. 

The second lesson comes from the resulting equation for the 
interface position. This equation takes the form of a weighted 
spatial average of the difference between the two transcrip- 
tion profiles (Eq. IS 14-b . In other words, in order to deter- 
mine the position of the interface, one needs to sum over 
many nuclei. In the sum, each nucleus contributes the dif- 
ference between the local transcription rates of the two RNA 
species. Clearly, these rates may be influenced by many fac- 
tors, and in a description which is somewhat closer to real- 
ity one would expect this difference to be fluctuating around 
a^(x) — a m (x). However, the interface position is a sum of 
these fluctuating objects, and one might hope that the sum of 
these fluctuations — which are uncorrected — would be close 
to zero. In this case, the interface position would be robust 
to fluctuations of this type. Indeed, a stochastic simulation of 
the model shows no change in the interface position (or struc- 
ture), as compared with the deterministic model discussed so 
far (see Supporting Information for details of our simulations). 

In passing let us note that the conditions required so far — 
namely, strong interaction between the miRNA and its target 
and small £ — may be reached by making the mRNA com- 
pletely stable (/3 m — ► 0). However, our analysis shows that in 
this case the system would never relax to a steady state, since 
target mRNAs would accumulate at the left end of the tissue 
without limit. Our analysis here is therefore only applicable 
if the mRNA molecules undergo autonomous degradation, in 
addition to the miRNA-dependent degradation. 



microRNAs can define sharp stripes of gene expression. 

Using the insight gained in the previous section we briefly 
show how a stripe is formed when the miRNA transcrip- 
tion profile a fl (x) is similar to a m (x) but displaced from 
it (Fig. EJ\). Suppose, for example, that the synthesis of an 
miRNA and its target are activated by the same transcription 
factor. In any given nucleus, the two promoters experience 
the same concentration of this transcription factor. However, 
they need not react in the same way: if the binding affinity of 
one promoter is stronger than that of the other, there will be 
intermediate concentrations of the transcription factor where 
the first promoter will be activated while the other will not. 
Such a scenario is depicted in Fig. 2A, where a common tran- 
scription factor, which exhibits a spatial gradient, activates the 
target gene as well as the miRNA gene. In this case, the target 
promoter has higher affinity to the transcription factor than the 
miRNA promoter. Thus, some cells in the middle of the de- 
veloping tissue express the target mRNA but not the miRNA. 

Unlike the case studied in the previous section, where the 
transcription profiles crossed at one point, here the transcrip- 
tion profiles cross at two points. Let us retrace our steps in 



the previous section, by first considering the case of no diffu- 
sion. For low values of the interaction rate k, the miRNA and 
mRNA profiles are qualitatively similar to their transcription 
profiles. As k is increased however, miRNA deplete mRNA 
levels at any position where a M > a m and thus confine 
mRNA expression to a stripe between the two crossing-points 
of the transcription profiles (green curve in Fig. [2)3). Can dif- 
fusion make this profile sharper, as in the previous case? In- 
deed, diffusing miRNAs which survive annihilation on the left 
and right diffuse into the interval between the two crossing 
points, and establish sharp interfaces in the mRNA concentra- 
tion profile. The resulting stripe resides within this interval, 
but is narrower (blue curve in Fig.[2]B). It is therefore impor- 
tant that parameters would allow for sharp interfaces, without 
making the stripe too narrow (or even disappear). Therefore, 
to sustain a well-defined stripe of gene expression, I, which 
determines the interface width, must be much smaller than the 
distance between the two crossing points of the transcription 
profiles. 

One can use the same analytic method mentioned earlier 
to calculate the new positions of the stripe boundaries (see 
Supporting Information). This exemplifies how the method 
can be used to analyze geometries of increasing complexity. 



Standard experimental methods may be used to probe the 
sharpening mechanism. 

The sharp interface that we predict can be detected di- 
rectly in an imaging experiment, provided the light inten- 
sity varies linearly with mRNA concentration, and that the 
spatial resolution is high enough. However an experimental 
readout is often not so faithful to the underlying concentra- 
tion profile. In the worst case, the apparatus' readout is bi- 
nary: concentrations below an apparatus-dependent threshold 
are not detected, while concentrations larger than this illicit 
a concentration-independent fluorescence intensity. In such 
scenarios, it is impossible to tell a smooth and sharp concen- 
tration profile apart as both yield a sharp interface in the binary 
readout (Figs. (3j\-C). In contrast, low spatial resolution may 
make a sharp interface appear smoother then it really is. 

Fortunately, our model of miRNA-mediated morphogenic 
regulation possesses another signature that is visible at such 
coarse experimental resolution. To detect this signature, one 
needs to overexpress the miRNA in a small patch of cells 
(hereafter denoted the 'patch'). Our model then predicts that 
this patch has a qualitatively different effect depending on 
which side of the interface it occurs. We confine ourselves 
in this section to the case of opposing transcription profiles of 
the miRNA and its target.. 

The technique one uses to generate the patch may differ 
according to the stage of development under consideration. In 
the early blastoderm stages of Drosophila development, for 
example, a Gal4 driver may be used to drive expression of the 
miRNA in those cells where an endogenous gene is expressed 
1 4411 . Many endogenous genes are expressed in stripes along 
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FIG. 2: Stripe of gene expression defined by microRNA interaction. A. Transcription profiles of a miRNA (red) and its target (green), favorable 
for stripe generation. See Supporting Information for parameter values. B. Expression of the target is restricted to the middle of the tissue even 
when the miRNA cannot diffuse (D = 0, green). Mobility of the miRNA makes target expression more focused (D — 0.001, blue): Target 
expression in the stripe is enhanced, and the stripe boundaries become sharp. 



the anterior-posterior axis during these stages and some have 
dedicated enhancers for single stripes Ir4l l45ll . As an example, 
the yeast FLP-FRT recombination system has been used to 
misexpress the gap gene knirps in a stripe by placing it under 
the control of the eve stripe 2 enhancer [4g| . 

In later stages of Drosophila development, e.g. imaginal 
discs, one technique is the random generation of a mosaic of 
mutant clones (patches) by mitotic recombination 144 14711 . 
The patches are generated at a low rate and one then screens 
for those embryos containing a single patch. 

We model the localized overexpression of miRNA by an 
effective local increase of the transcription rate by an amount 
a c = 5. This increase in transcription rate occurs in a small 
number of cells, which in our model is about 5% of the tissue 
length. More specifically, we choose a position x c for the 
center of the patch and, for every point x that resides within 
a distance w/2 of x c , we change the transcription rate from 
a^{x) to a^{x) + a c . Here w is the 'width' of the patch, 
which takes a value w ~ 5% of the tissue length. 

Consider positioning the patch first in the miRNA-rich re- 
gion of the developing tissue, Fig. 0D. One sees that, even 
if positioned at a distance from the expression domain of the 
target, the effect of the additional miRNA is to push the inter- 
face toward the left. The localized patch of cells therefore has 
a nonlocal effect. 

Intuitively, one identifies three contributors to this effect. 
First, of course, is the mobility of the miRNA. Diffusing miR- 
NAs are the carriers of information from the patch to the other 
cells. However, these miRNAs need to be able to travel ap- 
preciable distance before they are consumed, and thus the sec- 
ond factor is the fact that the patch was introduced in a region 
where miRNAs only experience autonomous degradation, and 
can thus travel a distance A (which we have assumed to be of 
the same order as the tissue length). Finally, the spreading 
miRNAs become localized again near the sharp interface with 
the mRNA-rich region. This is simply due to the fact that 
miRNA cannot penetrate this region over a distance that is 
larger than I (which we have assumed to be smaller than the 
width of a single cell). This intuitive picture is quantified in 
the Supporting Information. 

This experiment should be contrasted with one in which the 



overexpressing patch is positioned in the mRNA-rich region, 
as shown in Fig. [3}i. Such ectopic miRNA expression has a 
local effect, with excess miRNA creating a hole in the mRNA 
expression domain. The hole edges constitute two additional 
interfaces in the system, the sharpness of each again deter- 
mined by I. 

One can go further and make a quantitative prediction, re- 
lating the number of patches in a mosaic of patches with the 
lateral shift in the interface position. To first approximation, 
one needs to count the number of patches in the miRNA-rich 
region, and disregard completely the patches in the mRNA- 
rich region. The displacement of the interface position is then 
linearly proportional to this number; see Supporting Informa- 
tion for details. Simulated experimental results that would 
verify this prediction are shown in Fig.@] 

The distinct nonlocal effect described above does not occur 
when the miRNA are unable to move between cells. Also, 
we have checked (for the parameters used in this study) that 
it does not occur when the miRNA acts purely catalytically. 
Rather, both miRNA mobility and a strong interaction be- 
tween miRNA and target are required. The presence or ab- 
sence of the nonlocal effect would therefore confirm or falsify 
the hypotheses that miRNA are mobile and that they interact 
stoichiometrically with mRNA while in this mobile state. 



Irreversible localization of miRNA 

Finally, we consider the case where the miRNA transcrip- 
tion profile overlaps that of its target. This case is somewhat 
puzzling, as one would hardly see the point of going through 
this futile cycle. In particular, in the framework we have 
presented so far, such a coincidence of transcription profiles 
would reduce steady-state mRNA (and miRNA) levels to neg- 
ligible amounts in all cells. It is therefore not surprising, that 
this scenario is very atypical IT30ll . Where it does occur, e,g, 
possibly in the Hox gene network lf26ll . one explanation may 
be that this is only a transient state. In such a scenario, for 
example, the miRNAs keep the target mRNA from being ex- 
pressed, until such time when their translation is required. At 
that time, transcription of miRNA may be stopped, or alter- 
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FIG. 3: Indirect test of the sharpening mechanism. A. Concentration profiles without diffusion of miRNA (green) and with diffusion (blue). 
The black line denotes an apparatus-dependent threshold concentration. B, C. The readout of from an apparatus which amplifies any signal 
above its detection threshold does not distinguish between a smooth drop in mRNA levels (as occurs in B) from a sharp drop (the underlying 
profile in C). D. A patch (red) in the miRNA-rich region (posterior to the interface) pushes the interface to the left in a threshold assay. E. A 
patch in the mRNA-rich region (anterior to the interface) leaves the interface unaffected. See Supporting Information for parameter values. 




1 2 3 4 5 6 7 

number of patches posterior to interface 

FIG. 4: Quantitative test of the sharpening mechanism. Simulated 
experiment in which the interface position xt is measured in an en- 
semble of embryos having various numbers of patches. The theory 
predicts that the interface displacement should be linearly propor- 
tional to the number of patches. 



natively their interaction with mRNA may be inhibited. In- 
duction of the target would then proceed rapidly. An example 
where miRNA expression is regulated during the course of 
cellular development is found in T-lymphocytes 14811 : and reg- 
ulation of miRNA-target interaction has been demonstrated 
in neurons of Drosophila Ir49ll . Nevertheless, we consider 
here the case where the overlapping transcription profiles do 
represent — at least for some period of time — a steady-state, 
and show that with a slight change in the dynamics, this lay- 
out can still spatially separate the tissue into a region of high 
mRNA concentration and a region where target mRNAs do 
not accumulate. 

The key ingredient required to rescue mRNA levels from 
global obliteration is the possibility that only those miRNA 
localized to a cell interact with their target. Moreover, miRNA 
localization needs to be irreversible; the case of rapid re- 
versible localization reduces to the model studied in previous 
sections. 

The separation of transport and interaction may occur in 



a number of ways. miRNA may be transported from cell to 
cell in a vessel that needs to be dismantled before miRNA 
can interact with mRNA. Alternatively, the interacting form 
of the RISC complex may be immobilized. Another possi- 
bility is that interaction between miRNA and mRNA is lim- 
ited to specific loci in the cell, from which miRNA cannot 
escape. These scenarios — and others — are summarized by a 
model where miRNA acquire a mobile state immediately after 
synthesis, and can switch into an interacting localized state at 
later times 



(3) 



This model, a straightforward generalization of (Q~|i, is given in 
Eq. lS17l of the Supporting Information, where we also solve it 
exactly. 

Diffusion of the mobile miRNAs now occurs independently 
of the target, since mobile miRNAs do not interact with the 
mRNA species. Accordingly, a steady-state mobile miRNA 
profile is formed. This profile is smoother than the miRNA 
transcription profile (since this profile is formed by diffusion 
with no interaction). In particular, it is lower in those regions 
where the transcription profile a was higher than average and 
higher in those where a was lower than average; See Fig{5j\. 

The steady-state mobile mRNA profile is the source for lo- 
calized miRNA, which are thus generated in a spatially de- 
pendent manner. Their generation rate follows this profile, up 
to a factor q, the localization rate. This generation rate pro- 
file replaces the miRNA transcription profile of the previous 
sections. Thus, the role of the mobile miRNA is simply to 
set this profile. Left only with the localized miRNA and their 
mRNA targets, diffusion should not be considered any fur- 
ther. This problem has already been discussed in the context 
of opposing transcription profile, and from that discussion we 
realize that mRNA concentration still vanishes at any point 
where the miRNA generation rate profile exceeds that of the 
mRNA. In contrast, mRNA level is proportional to the differ- 
ence between the two profiles wherever the mRNA transcrip- 
tion rate is larger. Fig. 5B shows how the tissue is separated 
into two regions. We have shown that this scenario boils down 
to a diffusion-free problem, so the mechanism which could 
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FIG. 5: Irreversible localization rescues the target from annihilation 
throughout the tissue and sharpens its expression pattern. A. Steady 
state concentration profile of mobile miRNA (black) is proportional 
to the steady-state profile of the mobile miRNA, is a broadened form 
of the transcription profile a(x) (green). In this figure q = f3 s . B. 
As a result, mRNA (green) accumulates only at the left side of the 
tissue. The mRNA profile decays smoothly, and does not exhibit a 
sharp interface. 



sharpen the interface between the two regions is absent. As 
expected, the interface in Fig. 5B is smooth. 

The value acquired by the localization rate q is important 
for the success of the scheme described above. If this rate is 
too high, the resulting mobile-miRNA profile - and therefore, 
the localized miRNA generation profile - would coincide with 
the transcription profile, ultimately resulting in mRNA deple- 
tion in all cells. On the other hand, if the localization rate q is 
too low, the generation rate of mobile miRNA may be lower 
than the mRNA synthesis rate in all cells, and have very lit- 
tle effect on their accumulation. Taken together, we find that 
for 'optimal' performance, q should be set around the miRNA 
degradation rate, /3 M . However, our analysis shows, that no 
fine-tuning is required - in fact, setting q within an order of 
magnitude from /3 M would still yield two distinguishable re- 
gions along the developing tissue (Fig.[S4]and text in the Sup- 
porting Information). 



DISCUSSION 

In this study we have analyzed a model where miRNAs 
sharpen target gene expression patterns by generating an in- 
terface between high and low target expression. Such an in- 
terface can be directly detected provided the signal is not am- 
plified to saturation and provided the threshold for detection 
is sufficiently low. 

A necessary condition for the interface to occur is that the 
miRNA and target are co-degraded; a mode of miRNA-target 
interaction where miRNAs only promotes mRNA degrada- 



tion without being affected is insufficient. One can in prin- 
ciple test the existence of coupled miRNA-mRNA degrada- 
tion by inhibiting the transcription of mRNA and monitoring 
its decay rate, which in this case would be time-dependent, 
Pcs(t) = flm + ks(t). We note also that a stoichiometric 
interaction may complicate the interpretation of sensor trans- 
gene experiments M26I1 . as the transgene would then sequester 
miRNA and thereby alter the original expression patterns. 

The interface between low and high mRNA levels is char- 
acterized by low copy numbers of both RNA species. In such 
cases fluctuations in the copy number of either species may 
have macroscopic effects. For example, a small RNA-target 
pair in bacteria show enhanced fluctuations when their tran- 
scription rates become comparable |5(i 51]. These fluctua- 
tions can in turn give rise to noise-induced bistability, which 
manifests itself experimentally as diversity in a population of 
cells 115 111 . We performed Monte-Carlo simulations of our 
model, but found that fluctuations have no macroscopic ef- 
fect, even near the transition point where copy numbers of 
both species are low. This in-built robustness to fluctuations 
arises because the interface position is determined by an inte- 
grated transcriptional flux which averages out individual cel- 
lular fluxes. Thus spatial averaging results in high spatial pre- 
cision without smearing the interface. 

The model studied here is defined in one spatial dimension, 
which we identified as the axis x perpendicular to the inter- 
face between domains of gene expression. A two-dimensional 
extension of the model where the transcription profiles are 
independent of the y axis and diffusion occurs equally in 
both directions leaves the solution independent of the y axis. 
However, we have further generalized the model to allow the 
mRNA transcription profile to have a slight dependence on 
y (see Supporting Information). We find that the results de- 
scribed above hold also in this two-dimensional case. Further- 
more, we find that fluctuations in the interface position along 
the y-axis are considerably reduced by the miRNA regulation. 
A detailed analysis is deferred to a future study. 

Strong cooperative activation, as often occurs in morpho- 
genetic regulation at the transcriptional level (e.g. Bicoid has 
about five binding sites in target promoters of Drosophila), 
would seem to make pattern formation by morphogens inher- 
ently susceptible to temperature variations J52, 53]. Neverthe- 
less, embryonic patterning appears to be quite robust to tem- 
perature variations, as has been documented for Hunchback 



1 52|] and for Eve [54] in Drosophila. The only cooperative re- 
action required in the model presented in this work is coupled 
degradation of miRNA and mRNA, suggesting the possibil- 
ity that miRNAs filter fluctuations arising from temperature 
variations. 

Candidate systems in which to test the ideas put forth in 
this study include the establishment of dorsoventral (adax- 
ial/abaxial) leaf polarity in plants, as well as the segmentation 
of the early Drosophila embryo. We now discuss these two 
systems in some detail. 

Leaf polarity in plants is established shortly after the emer- 
gence of the leaf primordium from the meristem. Specifi- 
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FIG. 6: Model for leaf polarity in maize. A. In this model, TAS3 
is processed in the meristem into several ta-siRNAs. Two of those, 
ta-siR2141/2, migrate into the adaxial side of the primordial leaf, 
where it inhibits miR166 (either directly or via ARF3). miR166 goes 
on to repress expression of HD-ZIPIII genes, which define adaxial 
fate. B. Transcription profiles. Transcription of ta-siR2141/2 (red) 
occurs in the shoot meristem, next to the adaxial side of the devel- 
oping leaf. The target, either the miR166 or arf3, is transcribed uni- 
formly throughout the leaf (green). C. Steady-state concentrations of 
ta-siR2141/2 (red) and its target (green). Diffusion of ta-siR2141/2 
into the leaf restricts the expression of the target to the abaxial side, 
and - assuming a non-catalytic interaction - creates a sharp interface 
between the abaxial and adaxial domains. 



cation of leaf polarity depends on the Sussex signal 15511 . a 
meristem-borne signal that specifies adaxial cell fates. Mem- 
bers of class III of the homeodomain-leucine zipper (HD- 



ZIPIII) proteins specify adaxial fate 111 ll 15611 . In Arabidop 



sis and in maize, the polar expression pattern of these genes 
results from their inhibition by two miRNAs, miR165/166, 



which exhibit a complementary expression pattern III U 11211 . 
Recently it has been shown that in maize, restriction of 
miR 165/1 66 to the abaxial side of the developing leaf depends 
on the polarized expression of LBL-1, a protein involved 
in the biosynthesis of trans-acting RNAs, ta-siR2141/2 1 1311 . 



Possible targets of ta-siR2141/2 include members of the arf3 
gene family (a transcription factor that is ex pres sed abaxially) 
as well as members of the miR 166 family 111 311 . While miR- 
NAs in plants are thought to act mainly cell-autonomously 
|57|, DCL4-dependent siRNAs, such as ta-siR2141/2, may 



exhibit cell-to-cell movement 15811 . 

The following model is consistent with these data (Fig.|6j\). 
The RNA transcript TAS3 is cleaved to produce ta-siR2141/2 
in the meristem. These small RNAs then propagate (diffuse) 
into the adaxial side of the leaf, inhibiting the expression of 
miR166 either directly or through the ARF3 transcription fac- 
tor. The target (either miR166 or ARF3) is transcribed uni- 
formly throughout the leaf, and is localized to the cell where 
it is synthesized. If one further assumes that the interaction 
between ta-siR2141/2 and its target is non-catalytic, then this 
model belongs to the class of models studied in this work, 
and can therefore exhibit a sharp interface between the abax- 
ial domain of high target expression and the adaxial domain 
of no expression; see Fig. [6] In agreement with this model is 
the low abundance of ta-siR2141/2 in Arabidopsis |59[ 60ll . 
despite their distinct phenotypic role. 

Early embryonic development in Drosophila proceeds via 
a cascade of gene activities that progressively refine expres- 
sion patterns along the anterior-posterior axis of the embryo. 
A recent study of the expression patterns of nascent miRNA 
transcripts suggests that a number of miRNAs may play a role 
in this process. The miRNAs miR-309clus, miR-10 and iab- 
4 (which all reside between annotated mRNA genes on the 
genome), and miR-11, miR-274 and miR-281clus (which all 
reside within introns of annotated genes) are all expressed in 
a graded fashion along the anterior-posterior axis of the blas- 
toderm embryo ffl l6lll . 

The complementary transcription profiles of iab-4 and its 
target Ubx at stage 5 of development make this miRNA-target 
system a candidate for the sharpening mechanism proposed 
in this study. The early Ubx transcript pattern is broadly dis- 
tributed over the posterior half of the embryo, becoming lo- 
calized to a stripe at the centre of the embryo by the com- 
pletion of cell formation |2^, 62], probably as a result of tran- 
scriptional repression by Hunchback in the anterior and poste- 
rior regions of the embryo |63l] . The nascent transcript profile 
of its regulator, iab-4, is broadly distributed posterior to this 
stripe 12911 . 

It may be the case that iab-4 is also expressed before cell 
formation and that the absence of cell walls makes iab-4 mo- 
bile. The much larger Ubx mRNA, on the other hand, may 
be effectively stationary on the timescales of interest |7|]. Fur- 
thermore, the transcription profiles of iab-4 and Ubx at stage 5 
do not seem to overlap l29ll . suggesting that iab-4 intercellular 
mobility may be crucial to allow it to interact with Ubx at this 
stage of development. Assuming then that only the miRNA 
iab-4 is mobile, the complementary expression patterns of iab- 
4 and its target Ubx measured in Ref. I29I1 is consistent with 
our model of miRNA-induced sharpening. Sample profiles 
predicted by the model are shown in Fig. [7] 

A possible difficulty with regard to applying our model to 
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FIG. 7: Model for sharpening of Ubx expression boundaries by iab- 
4. A. The assumed transcription profiles of Ubx (green) and iab-4 
(red).B. In the absence of iab-4 (or when the interaction between iab- 
4 and Ubx is inhibited), the expression pattern of Ubx is broad and 
smeared. C. When iab-4 and Ubx interact, the expression pattern of 
Ubx becomes narrow and sharp. The diffusion constant is D = 0.1 
and the co-degradation rate is A: = 1. The embryo-like shape is for 
illustrative purpose only. 



Ubx/iab-4 is that the system may not have reached steady state 
before stage 6 when cells begin to migrate. In particular, no 
Ubx protein was detected at stage 5, possibly because of the 
time needed to transcribe the large Ubx locus l29tl . 

Like iab-4, the miRNA miR-10 is also expressed at stage 
5 in a broad posterior region along the anterior-posterior axis 
fl 14f | . The homeotic gene Scr is a predicted target of miR- 
10 |64] and is also expressed in the blastoderm at stage 5 16511 . 
The miR- 1 site in the Scr 3 ' UTR is likely to be functional be- 
cause the pairing is well conserved in all drosophilid genomes 
and because the miRNA site is conserved in the Scr genes 
in mosquito, the flour beetle and the silk moth l20Tl . Unlike 
Ubx, the protein of Scr is detected at this stage of develop- 
ment in a stripe of ectodermal cells about 4 cells wide in the 
parasegment-2 region, though it may not be functional at this 
time as the protein (a transcription factor) was not localized to 
the nucleus 16511 . This spatial expression pattern is proximal 
to the anterior limit of miR-10 expression 



65H . Hence the 



interaction of miR-10 with Scr at stage 5 of Drosophila devel- 
opment is also a candidate for the sharpening mechanism. 

The sharpening mechanism is most effective when the spa- 
tial transcription profiles of miRNA and target are regulated 
in such a way as to be mutually exclusive. The genomic loca- 
tions of the miRNAs iab-4 and miR-10 are proximal to their 
targets, which is certainly consistent with the possibility of 
coordinated regulation 16611 . 
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SUPPORTING INFORMATION 



Parameter Values 



1. The transcription 



Throughout, j3 m = [3^ L = 0.01, k 
profiles of Fig.Q]are 

a m (x) = 0.5A m [tanh((a; tsa! -x)/ \ tsx ) + 1] (SI) 
a^(x) = 0.5^[tanh((a; - x tsx )/X tsx ) + 1], 

where A m = 2, = 1, xt sx = 0.5 and A tsx = 0.2. In 
the stripe geometry (Fig. [2]), the transcription profile for m is 
as above, with xt sx — 0.3. The transcription profile of the 
miRNA, is given by 

a M (x) = A^o + 0.5A AI [tanh((a; tSK - x)/\ tsx ) + 1], (S2) 



with A, 



2, A, 



iy, — ~, ^i M o = 0.6 and x tS2; = 0.7. To test our model 
for irreversible localization we use the same transcription pro- 
files as in ( fsTb . with A m = A^ = 1. 

In the discussion we outline possible applications in two 
systems: leaf polarity in maize and segmentation in the early 
Drosophila embryo. Here we did not aim to estimate parame- 
ters from experimental data (which, in most cases, is not quan- 
titative enough to allow for parameter inference). Instead, 
parameters were chosen arbitrarily to allow clear demonstra- 
tion of possible results. In the case of leaf polarity (Fig. [6} 
we choose a m (x) = A m ,a s (x) = A ll 6{x - x tsx ) with 
A m = IjAfi = 50 and xt sx = 0.99. Here Q(x) is the unit 
step function. In the Drosophila embryo (Fig. |7J, the tran- 
scription profile of the iab-4 miRNA is the same as in (IS 1 
while the Ubx mRNA transcription profile is given by 



oi m {x) = 0.5A m [ta,nh(\xt sx - x\/X tsx ) + 1] 



(S3) 



To study a two-dimensional extension of this model, we let 
the mRNA transcription profile have a slight dependence on 
the y-direction, e.g. by setting xt sx in ot m (x) (Eq. IS It . to 
x tsx = sin(47ry)/10. 

Strong-interaction limit 

Here we define, for the case of no diffusion, what we mean 
by the strong-interaction limit, but the limit is equally valid 
when diffusion is present. We make the critical assumption 
that at any cellular position the decay of the minority species 
is dominated by coupled degradation; by rescaling Eq. [2] of 
the main text, this condition amounts to 

•fccUsVftn ! K i MZM»l, forallx , (S4) 
Pn P™ J 

In the absence of diffusion, Eq. [2] boils down to a quadratic 
equation for the steady-state mRNA level m, and one can im- 
mediately write the solution 



m = f3j 



,(1 + e) + ^/[a m + a M (l + e)] 



(S5) 



where e = ^^/(^a^), which is assumed in flS4l to be 
small. To zeroth order in e, this expression simplifies to 



[ot m - a M ] + 
m« , 

where [x]+ = max{0, x}, as depicted in Fig.[Tfc. 



(S6) 



Analytical Approximation 

To understand the origin of the length scales A and I, and 
their relation to the tissue length L when the interface is sharp, 
consider first the region of space where miRNA are in the 
majority. In this region, where kfi (3 m , we neglect the 
autonomous-degradation term in d2al ). yielding a m = kmfi 
and thus 



= dtp — a r 



PuV + D(i" 



(S7) 



Hence miRNA are produced at an effective rate — a m and 
diffuse over distances of order 



A = J — 



(S8) 



which, as we have argued in the main text, should be com- 
parable to the tissue length, A ~ L. On the other hand, in 
the mRNA-rich zone km ^> (3^ and so the only length scale 
available to the miRNA is 



D 



fccC/pV) 



(S9) 



where a* n is a typical value of a m (x) in the mRNA-rich zone. 
For the mRNA this is the only length scale that competes with 
the spatial layout provided by the transcription profile and so 
it must determine the interface width up to a constant prefac- 
tor p. Hence the second condition for a sharp interface is that 
pi <C L. In agreement with the expression in Eq. [S9] de- 
rived on hueristic grounds, our numerical solutions show that 
the interface becomes broader when the co-degradation rate k 
is decreased (Fig.QJj) or when the diffusion constant D is in- 
creased (Fig.Q]H). We note that the limit p£ <C Lis equivalent 
to 



ka* m /p m > D/L 2 



(S10) 



In other words, a sharp interface arises when the co- 
degradation rate of miRNA and target dominates the rate of 
diffusion over macroscopic distances. Since we are neglect- 
ing diffusion of miRNA in the mRNA-rich region, the mRNA 
profile there is given again by (|S6l l. 

' Microscopically, miRNA in the miRNA-rich zone dif- 
fuse in a landscape dominated by autonomous degradation, 



2 



leading to the decay length A in (IS8b . Upon entering the 
mRNA-rich region, co-degradation suddenly overwhelms au- 
tonomous degradation of miRNA (km 3> /3„), and the 
miRNA are faced with an effective absorbing boundary. We 
therefore expect the miRNA concentration to vanish as one 
approaches the interface from the right. In addition, our pic- 
ture asserts that the miRNA concentration is vanishingly small 
everywhere on the left of the interface. Taken together, these 
two properties impose zero miRNA concentration and zero 
miRNA diffusive flux at the interface. These two boundary 
conditions on the miRNA dynamics at the interface between 
mRNA-rich and miRNA-rich regions, together with the zero- 
flux condition at x = L, allow us to determine the position 
Xt of the interface. Furthermore, the interface must lie in the 
region defined by a m > because co-degradation can dom- 
inate autonomous degradation of miRNA only if there is a 
reservoir of mRNA to co-degrade with. 

Armed with this insight into the miRNA profile we 
may now solve dS71 i, subject to the boundary conditions 
—Dfi'(xt) = —Dfx'(L) = and /i(x t ) = 0, in terms of a 
Green's function. Making use of the zero-flux boundary con- 
ditions, the Green's function of (IS7b is 



g(x,s) = 



G(x, s) if x < s 
G(s,x) if x > s ' 



(SH) 




FIG. SI: A comparison of the exact numerical concentration profiles 
with the analytical approximation (grey lines) described in the text. 
Note that the interface is always sharp in the analytic approximation. 
The diffusion constant here is D — 0.01 and the miRNA-mRNA 
interaction parameter is = 1. Autonomous degradation rates are 
fa = 0.01. 




where 



40 60 

cell index 



100 



XG{x,a) 



cosh (2=2*) cosh (±=i) 



sinh 



' L—Xf 



(S12) 



The miRNA profile is then a weighted spatial average of the 
net transcriptional flux of miRNA to the right of the interface 



FIG. S2: Temporal averages (squares) of a stochastic simulation 
of the model, compared with the corresponding mean-field solution 
(solid lines), in a developing field comprising 100 cells. The diffu- 
sion constant is D = 1000 and the miRNA-mRNA interaction pa- 
rameter is k = 100. The autonmous degradation rates are pi = 0.1. 



[a M (s) - ce m (s)]g(x, s) ds. 



(SB) 



Employing the zero-concentration boundary condition 
n(xt) = 0, we arrive at the following implicit equation for xt 



Stochastic Simulations 



[a M (s) - a m (s)] g(x t , s) ds = 0. 



(S14) 



Solving for x t requires knowledge only of the transcription 
profiles. Moreover, one immediately sees that x t can tolerate 
fluctuations in the transcription profiles which preserve the in- 
tegral. This should be contrasted with the non-diffusive case 
in which xt reduces to the crossing point of the transcription 
profiles, a m (xt) = atp(xt), which is less robust to small- 
number fluctuations. 

In Fig.lSllwe compare the analytical expressions in Eqs. lS61 
IS13l and lS14l with the exact numerical solution of 10. As ex- 
pected the agreement is good because co-degradation domi- 
nates both autonomous degradation (Eq. IS4b and loss due to 
diffusion (Eq.HTO). 



We used the Gillespie algorithm to stochastically simulate 
the reaction and diffusion events on a one-dimensional grid 
of cells fl67ll . In this algorithm the next event, as well as the 
time to the next event, are chosen randomly. A simulation 
using 100 cells (the approximate anterior-posterior length of 
the Drosophila embryo during cycle 14) is compared with the 
corresponding deterministic solution in Fig.[S2] Surprisingly, 
the deterministic solution is a good approximation for mRNA 
anterior abundances as low as 20 molecules per cell. As ex- 
pected, in cases where the predicted interface is of the order of 
a single cell we find that the solution to our model underesti- 
mates the width of the interface. This, for example, is the case 
when the developing tissue becomes as small as 10 cells — the 
approximate size of the leaf-organ primordium during plant 
development. 
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FIG. S3: A comparison of the exact numerical concentration profiles 
with the analytical approximation (grey lines) described in the text 
when the transcription profiles are aligned so as to produce a stripe. 
See caption of Fig. [ST1 for parameter values. 

Stripe Boundaries 

We denote the left-most boundary of the stripe by Xti and 
the right-most by Xti- The location of these interfaces is de- 
termined analytically by solving (IS7b with zero-flux bound- 
ary conditions in the interval [0, Xti) (and then enforcing 
n{xti) — 0) and in the interval [xt2, L] (and then enforcing 
n{xt2) = 0). In the region between the interfaces the mRNA 
profile is approximately given by m = \a m — a^]/ ' fi m ; in 
the portion of the developing tissue complementary to this the 
mRNA profile is negligible. The analytic profiles for /i and m 
are compared with the exact numerical solutions in Fig. IS3I 



Experimental Prediction 

The nonlocal effect predicted when the miRNA is overex- 
pressed in a small patch of cells can be understood with the 
aid of ( |S7t . First, integrate this equation from the interface to 
L, both with (i t ) and without (x t ) the patch. Then, neglecting 
the change in the miRNA concentration induced by the patch, 
one can show that 




(S15) 



In other words miRNAs from the patch diffuse across the tis- 
sue towards the region [x t ,x t ] and annihilate all mRNA pro- 
duced there that would otherwise have maintained the inter- 
face at Xt. 

As mentioned in the main text, we can make a quantita- 
tive testable prediction when there are a number of indepen- 
dent patches. To see this recall that Eq. IS 151 is an integral 
relationship. Hence, in the corresponding equation for multi- 
ple patches, the left-hand side is proportional to the number 
of patches, provided we make the reasonable assumption that 
the patches are uniform in size and transcription rate. Simi- 
larly, the right-hand side is proportional to the interface shift, 
xt — it, provided the mRNA transcription profile is suffi- 
ciently flat in the interval [i t , xt]. Hence the interface position 
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FIG. S4: The ratio between mRNA concentration at the left (x — 0) 
and right (a; = 1) ends of the tissue, for the transcription profiles and 
parameters of Fig. [5] In the absence of miRNA regulation, this ratio 
is ~ 3.9 (horizontal dotted line). microRNA regulation is most 
effective in clearing residual mRNA from 'off' cells when q ~ /3 M 
(vertical dotted line). 

it decreases linearly with the number of patches. This predic- 
tion can readily be tested using an ensemble of embryos with 
varying numbers of patches. Simulated experimental results 
that would verify this prediction are shown in Fig. [4] in the 
main text. 

Irreversible Localization 

We employ this flavor of the model for the special case that 
the miRNA and mRNA transcription profiles coincide. Mo- 
bile miRNA can no longer interact with their targets; only 
when they flip irreversibly to the localized form 

fi -^-> fj, e . (S16) 

can they do so. This adds a new equation to our model which 
now takes the form 

= a — /3 m m — km(j,t (S17a) 
= q/j, — Ptfie — kmm (S17b) 
= a- Ppii-qii + Dn" . (S17c) 

Here fj,£ is the concentration of miRNA in the localized state, 
and q is an effective rate at which mobile miRNA switch into 
the interacting localized mode. The common transcription 
profiles of mobile miRNA and mRNA is denoted by a(x). 
The mobile-miRNA profile fi(x) is determined only by the 
third equation and, by analogy with Eq. IS131 is a weighted 
spatial average of a. As such the effective transcription pro- 
file of immobile miRNA q/i(x) is reduced on the left and in- 
creased on the right relative to that of the target a(x) The first 
two equations then reduce to the zero-dimensional case con- 
sidered earlier in which miRNA clear up mRNA to the right 
of the point where qfJ,(x) and a(x) cross. 

The value of q which is most favorable for the failsafe 
function is q w To understand why, it is helpful to 

consider two limits. In the first, localization dominates all 
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other loss processes for mobile miRNA, i.e. q 3> /3 M and 
q 3> D/L 2 . Hence all mobile miRNA immediately become 
immobile, the transcription profiles of mRNA and immobile 
miRNA coincide (a « q/i, Eq. IS 17ct , and all mRNA are im- 
mediately degraded in all cells. In the second limit, local- 
ization is negligible compared with independent degradation, 
q "C /3 M , and very few mobile miRNA become immobile. 
Consequently q[i <C a and there are practically no immo- 
bile miRNA available to affect the target. One measure of 



the efficacy of the failsafe mechanism (threshold response) is 
the ratio between mRNA concentration at the two ends of the 
tissue, m(0)/m(L). With the transcription profiles and pa- 
rameters of Fig. [5] this ratio would be ~ 3.9 in the absence 
of miRNA regulation, and ~ 257.4 in its presence. In Fig. IS4I 
we plot this ratio as a function of q, keeping all other parame- 
ters as in Fig. |5] One sees that the failsafe mechanism is most 
effective in the intermediate regime where q « (3^. 



